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RECOGNIZING AND DEVELOPING CREATIVE POTENTIAL* 


JOHN E. ARNOLD 
Professor of Mechanical Engineering, Massachusetts Institute of Technology 


OU CANNOT TELL, by examining the end product, 

whether it was obtained by means of creative activity or 
by means of more prosaic, routine methods, or even by pure 
chance. And, you cannot tell, by looking at the types of 
problems that men set for themselves, whether they are cre- 
ative or noncreative. The great scientists of the past and the 
present are great because they attacked their analytical prob- 
lem with great daring and imagination. ... Unfortunately, it 
is also possible to arrive at satisfactory answers to the multi- 
solutional, creative type of problem by extremely prosaic, 
nonimaginative, noncreative means... . 


Attributes of the Creative Individual: Now, what are 
some of the attributes that define or describe the creative in- 
dividual? In general, they fall into two categories—mental 
and emotional. The creative person is more “sensitive to 
problems” than the noncreative person. He can see need 
areas that require attention long before the less sensitive per- 
son. He not only has the ability, but he enjoys searching in 
the dark areas and leaves the work in the well-lighted places 
to his less imaginative brethren. He has developed a con- 
stantly questioning attitude, a healthy skepticism, rather than 
a complacent acceptance. The questions he asks and the 
problems he poses for himself are stated in such a way that 
they do not limit or confine his thinking or actions. Correct 
problem statement is basic to successful creative activity. In- 
correct problem statement can precondition the thinking and 
preclude desirable solutions. It also can have a tremendously 
inhibiting effect on three other mental attributes: fluency, 
flexibility and originality. 


Fluency of Thought: The creative person is more fluent 
in his thinking than the noncreative person. He has more 
ideas per unit time than others. Obviously, they are not all 
good ideas. But when he completes his list of alternatives, 
the probability is very high that there will be one or more 
ideas that are novel and that will solve his problem better. 
While it is very doubtful that quantity breeds quality, my 
experience has been that those who have lots of ideas, usually 
have the best ones. 


Flexibility of Thought: Flexibility of thinking reflects 
itself in the wide variety of approaches that the creative per- 
son chooses to investigate. The noncreative person’s past 
experience provides him with a comfortable little rut in 
which to operate, and he has great difficulty getting out of 
that groove. 

One of the worst examples of rigidity of thinking that I 
see as I go around and talk to and work with various indus- 
trial groups is the way in which they limit the design and 
research activity to only one or two of the four broad areas 
in which design activity can take place. Products or machines 
can be improved by increasing their function, by making the 
item do more than it did before, or by making it accomplish 


* From a pepe presented before the ASME Annual Meeting, Chi- 
cago, Nov., 1955, and appearing in Machine Design, May 3, 1956. 


its prime purpose by an entirely new approach. The per- 
formance level of a product can be raised by making the 
product longer lived, more reliable, more accurate, safer and 
more convenient to use, and easier to repair and maintain. 
The cost may be lowered by eliminating excess parts, substi- 
tution of cheaper materials, cheaper and more efficient manu- 
facturing methods, designing for convenience of assembly, 
designing to reduce hand labor or for complete automation. 
Last, products can be improved by designing for increased 
salability, that is, improving the appearance by making the 
product, its package and its point of sale more attractive, and 
by having a better appreciation of public wants. 


Originality: The creative person thinks of more unusual, 
more unique solutions to problems than the noncreative per- 
son. He apparently has the ability to penetrate very deeply 
to remote regions. He can break down the seemingly natural 
boundaries within his own mind, so that his experience, in- 
stead of being pigeon-holed, is completely and carefully 
cross referenced. The creative person, then, is more sensitive 
to problems, he is more fluent and flexible in his thinking, 
and he is more original in the solutions he arrives at. We 
might say, he is better able to question, observe, associate 
and predict. 


Emotional Attributes: The mental attributes listed are 
necessary, but not sufficient to the creative personality. 
Along with them there are a number of emotional attributes 
that many, I am sure, would list as primary rather than sec- 
ondary. First of all, a man must be motivated before he will 
begin to attempt to solve any kind of problem, analytical or 
creative. We have ample evidence that over-motivation is 
extremely disastrous, in that it inhibits creative activity, but 
I am afraid that more of us suffer from under-motivation 
rather than over-motivation. Passivity is probably the orig- 
inal sin, and initiative is what most of us lack, for once into 
a problem, most of us can generate enough interest to keep 
us going until the problem is solved. 

Willingness to take a chance, to gamble, may be another 
form of initiative. But it is an important aspect of the cre- 
ative person’s personality. The creative person has to be 
daring. He has to be a leader in his group for society, and 
he must constantly take calculated risks in his attempt to 
find better solutions to the problems that face mankind. He 
cannot stick to the safe, the tried and true, the prosaic ap- 
proaches, and he must pioneer in new areas in a very daring 
fashion. Creating, unfortunately, also involves destroying. 
The man who is seeking a new, better solution to an old 
problem is doing so, in part, because he wants to destroy a 
present, possibly adequate solution. 

Properly motivated and willing to take a chance, the cre- 
ative worker must, in addition, have self-confidence in his 
own ability to come up with a new and better solution. This 
is an extremely important emotional attribute and can only 
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CONSIDER A CAREER IN ENGINEERING TEACHING* 


HAROLD E. WESSMAN 


Dean of the College of Engineering 


A tremendous amount of 
publicity in recent months has 
focused attention on the short- 
age of engineers and on the 
competition by industry and by 
government agencies for young 
graduates. Every engineering 
college campus is overrun by 
trained recruiters coming to in- 
terview the seniors and the 
graduate students. Stories about 
outstanding students being 
sought by twenty or more in- 
dustrial concerns are becoming 
commonplace. Is it possible 
that in this intense competi- 
tion by industry and by gov- 
ernment, the opportunities for careers in teaching are being 
overlooked by our young engineers and scientists? 

Engineering colleges do not, and cannot, engage in the 
highly organized recruiting activities characteristic of large 
industrial organizations. Nevertheless, enrollments in engi- 
neering colleges are increasing rapidly and many new 
teachers are needed immediately. Many more will be 
needed in the years ahead. And so, to the young seniors 
and graduate students with outstanding records, I say, “Do 
not overlook the fine opportunities open to you in engineer- 
ing education. Go to your department head or dean and 
learn how stimulating it is to be an engineering teacher, to 
work with young men and to help them get the knowledge 
needed by a professional engineer.” 

Statistics published by the U. S. Office of Education show 
that 243,000 students were enrolled in our engineering 
colleges in the Fall term 1955. Analyses of population 
statistics and college trends have resulted in a prediction 
that this engineering enrollment will double by 1967. One 
might readily conclude that we will need twice as many 
teachers in 1967 as we have today. However, the increase 
in graduate study and research programs wherein the ratio 
of teachers to students is greater than it is for undergadu- 
ate programs means that we will need more than twice as 
many teachers. 

A recent survey made by the Relations with Industry 
Division of the American Society for Engineering Educa- 
tion reports that there are now approximately 8500 engi- 
neering teachers in the United States. However, some of 
these men retire each year because of age, and there are 
other losses due to health or illness. Consequently, about 
10,000 new engineering teachers will be needed in the next 
ten years. In fact, the ASEE survey just mentioned states 
that 1300 new teachers are needed for this academic year 
1956-57. Obviously, there is no lack of opportunities. 


H. E. Wessman 


Lists qualities of a good teacher 

What are the qualities of a good teacher? Of course, he 
must know his subject, know it well—so well that he can 
restate in his own words what others have written into the 
textbooks he may use in his classes. Moreover, he must 
know a great deal about supporting areas of knowledge. 
But do not conclude that a mastery of the knowledge in a 


* Reprinted from Industrial Science and Engineering, October, 
1956. 


particular field must be achieved before one enters teach- 
ing. Every teacher knows that he learned more about his 
subject during his first year of teaching than he did during 
all of his days as a college student. In fact, if he is a good 
teacher, he never stops learning. The man in industry, 
subject to the pressures of an industrial routine, does not 
usually have the same opportunities to achieve a mastery of 
knowledge as the man in a university environment. 

But there is even more to teaching than knowing one’s 
subject. A good teacher speaks clearly. His explanations 
are understandable. He shows interest and enthusiasm in 
his subject and he gets students interested in the subject. 
He takes an active, personal interest in the progress of his 
class. He is friendly and sympathetic. A good teacher has 
good manners, good morals, and a good sense of humor. 

What are the elements which make teaching such a 
stimulating and satisfying experience? Working with young 
men and watching their minds develop is one of the great 
sources of satisfaction. In fact, every teacher takes great 
pride in watching the students who have sat in his classes 
go out into the world and become successful men and 
citizens. Every teacher feels that he has contributed some- 
thing to that accomplishment. 


Teaching is stimulating 

Another inviting factor in teaching is the stimulation 
which comes from a thorough understanding of the lim- 
itations of existing knowledge and the acceptance of the 
challenge to discover new facts through research. There 
are many gaps and unexplored areas in the accumulation 
of scientific and technical knowledge. Those individuals 
in a university environment have unusual opportunities to 
explore the area of the unknown and to bring light and 
truth where ignorance prevailed before. It is difficult to 
describe the joy, the satisfaction, experienced by those who 
have discovered new principles, or new formulae, or new 
and better ways to solve problems, or new information 
about engineering materials and technical processes. 

Another element which makes college teaching so at- 
tractive is the association with kindred teachers and research 
workers in many diverse fields. The university is a com- 
munity of scholars engaged in creative work in teaching, 
writing, research, or in professional activity. It is difficult 
to find any other career with the same inspiring environ- 
ment. Moreover, the opportunity for a fine family social 
life in a university community is another good reason why 
many people find happiness and satisfaction in teaching 
careers. e 


Salaries on the upswing 

How about salaries in engineering teaching? They are 
not as high generally as salaries in industrial positions. In 
fact, the salary pattern in many schools is quite unsatisfac- 
tory at present. But steps are being taken to remedy this. 
Many college administrators now recognize that salaries in 
engineering teaching must be brought to the level of 
teaching salaries in medicine and law, if the quality of tech- 
nical education in the United States is to be maintained and 
improved—in short, if we are to survive in the technolog- 
ical race for world supremacy. Basic salaries in engineering 
teaching are on the upgrade. Moreover, extra income may 
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LEFT 


Fic. 1. Upstream view of 
model corresponding to sketch 
shown in Part I. 


BELOW 


Fic. 4. Using dye to trace 
flow distribution. 


RIGHT 


Fic. 2. Downstream view 

of model. Fish ladder at left; 

power-house collection system 

at right; forebay in upper right. 


ABOVE 


Fic. 3. Measuring water 
elevations throughout 
the model. 
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A Model Study of Some Fish Facilities 
for the Proposed Mayfield Dam 
Part I: Upstream Migrant System* 


E. P.. RICHEY 
Associate Professor of Civil Engineering 


The basic intent of the 
upstream system is to at- 
tract the migrants from the 
river into a collection sys- 
tem and provide a ladder 
for raising them to an ele- 
vation above the highest 
water level where they may 
be trapped for mechanical 
conveyance to upstream re- 
lease points. Reference is 
made to the January issue 
of The Trend for a sketch 
of the proposed dam and 
fish facilities. The barrier across the river will pre- 
vent the fish from bypassing the entrances to the 
collection system. These entrances are provided in 


E. P. Richey 


the powerhouse between the turbine draft tubes as | 


well as near the junction of the powerhouse and the 
barrier. The ladder begins at the end of the power- 
house and makes a 180° bend, finally to terminate in 
the traps at the base of the “upstream facility.” 

Figure 1 shows an over-all photograph of the 
model, built to a scale of 1:10, and was taken from 
‘the same angle as shown in the sketch; in the proto- 
type, the section in the foreground would be con- 
tained within the powerhouse. The entrances to 
the fish-collection channel are shown discharging 
into the tailrace. The main entrance and the barrier 
appear in the upper right corner, and the ladder 
portion occupies the background. A view of the 
model with the opposite perspective is shown in 
Fig. 2. The ladder ends near the manometer board 
in the center of the picture. (See also Fig. 3.) 

In the hydraulic sense, the upstream facility is a 
complicated open channel that is required to main- 
tain certain standards of flow throughout the oper- 
ating range of tailwater elevations. The amount of 


* The Upstream Migrant System, another part of the fish 
facilities for the proposed Mayfield Dam on the Cowlitz 
River, was studied in model scale at the Charles W. Harris 
Hydraulics Laboratory. This article is a sequel to Part I: 
“The Skimmer” that appeared in the January issue of The 
Trend. Both studies were carried out on models designed 
and built by the Major Projects Division of the Tacoma De- 
partment of Public Utilities. 
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H. H. CHENOWETH 
Assistant Professor of Civil Engineering 


water delivered must be 
varied with tailwater in 
order to keep the velocities 
in the transportation system 
reasonably constant. For 
instance, at high tailwater, 
the correct quantities would 
have to be _ introduced 
through diffusion gratings 
into the transportation sec- 
tion at seventeen different 
points. The water for oper- 
ating the upstream system is 
to be pumped from the river 
rather than drawn from the reservoir behind the dam. 
Thus, the amount of water necessary and the eleva- 
tion it must be lifted became important items in the 
design. The purpose of the model was to check the 
initial design layout and to observe flow conditions 
for different arrangements of the basic components 
of the system. 


The specified design flow rates must be discharged 
from the several entrances in order to attract the 
fish into the collection channel. When the fish are 
once in the channel, the depth and velocity of the 
water must be conducive to their upstream move- 
ment. In order to maintain the required velocities in 
the powerhouse collection channel with a minimum 
amount of water, various types of baffles were used 
in the powerhouse section. The model indicated the 
flow pattern (Fig. 4), making a comparison of the 
different baffles relatively easy. One of the series of 
baffle arrangements is shown in Fig. 5. 


H. H. Chenoweth 


The powerhouse collection channel and the main 
entrance bring the fish to the ladder. Each pool in 
the ladder series shown in Fig. 6 provides a one-foot 
rise in elevation. The fish may pass over the notched 
weir or through a port in the weir plate, not visible 
in the photograph. In order to provide additional 
water to maintain the proper velocity with increasing 
tailwater elevation, a supply channel is provided be- 
tween the two legs of the U-shaped ladder. The 
added water flows from the supply channel over the 
vertical plates, or “chimneys,” on the side of the 
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BELOW 


Fic. 6. One of the ladder series 
showing orifices. 


LEFT 
FiG. 7. Initial design of the ladder bend. 


Fic. 5. View of one type of baffle. 


RIGHT 
Fic. 8. Final design of the ladder bend. 
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channel, then downward and through a grating to 
discharge laterally into the ladder proper. In the 
initial design, the flow through the chimneys as well 
as through other parts of the system was not uni- 
form. A basic revision was made by inserting flow- 
metering orifices between the supply channels and 
the collection and transportation sections. The ori- 
fices were dimensioned to increase the head differen- 
tial needed to deliver the desired flow rates from a few 
hundredths of a foot to one foot. The resulting rates 
were thus made much less sensitive to minor fluctu- 
ations in surface elevations. The orifices between the 
supply channel and the ladder also are shown in 
Fig. 6. Similar orifices were located upstream from 
all other diffusion gratings. The orifices can be placed 
at different elevations to function quite automatically 
to meet the variable flow demands as the tailwater 
elevation changes. Their stabilizing effect on the 
delivery of water to the various parts of the facility, 
however, comes at the sacrifice of increasing the head 
against which the pumps must operate. 

The model did not disclose any new hydraulic con- 
cepts, nor was this its intent. It did provide a rapid 
way of showing how a wide variety of different con- 
trol features could be brought together which would 


be impractical to attempt by theory alone. A fish 
facility presents very important aspects of flow con- 
ditions that are of little consequence to the usual 
hydraulic engineering problems. One example of this 
is brought out in the two photographs of the bend 
in the ladder. In the initial design (Fig. 7) an up- 
welling occurred on the curved wall. So far as the 
pure hydraulics is concerned, this feature is not sig- 
nificant. However, such a current adversely attracts 
fish to delay their upstream migration, so it must be 
corrected. A reorientation of the transverse weir 
plate to the position shown in Fig. 8 minimized the 
undesirable upwelling. 

Another feature, undesirable from any viewpoint, 
would be the development of a surface wave train 
extending through the notched weir plates. Attempts 
to incite such a wave pattern in the ladder section 
were unsuccessful, leading to the conclusion that the 
types of weirs installed and the relative dimensions of 
the ladder would inhibit the formation of such a wave 
pattern. Thus, the model furnished a means of meet- 
ing the need for visual and qualitative information on 
the over-all performance as well as the usual quanti- 
tative measurements to compare the performance with 
the design standards. 


AEC Grant to Provide Nuclear Equipment 


HE COLLEGE OF ENGINEERING has been granted 
$109,500 by the Atomic Energy Commission for 
the purchase of nuclear laboratory equipment and 
materials. The grant will be used to equip the Re- 
actor Theory, Chemical Engineering, Nuclear Met- 
allurgy, and Radioactive Waste Disposal laboratories. 


The major unit in the Reactor Theory Laboratory 
will be a subcritical graphite exponential assembly 
similar to the one at the Argonne National Labor- 
atory. Parts for a reactor simulator, which will be 
built at the University of Washington, will also be 
obtained with the grant. 


The major part of the grant will be used for a 
number of small items, such as scalers, counters, and 
survey meters with which all the laboratories must 
be equipped. The Chemical Engineering Laboratory 
will also be equipped with a spectrophotometer, a 
scintillation spectrometer, a sample changer, and 
filter hoods. 


Equipment for the Nuclear Metallurgy Laboratory 
will be expanded to include a metallograph, a dil- 
atometer, a vacuum furnace, a heat-treat furnace, and 
an abrasive cutoff wheel. 
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A graduate program in Nuclear Engineering lead- 
ing to a Master of Science degree was initiated in the 
fall of 1956. It is a cooperative undertaking of the 
Departments of Chemical, Civil, Electrical, Mechan- 
ical, and Metallurgical Engineering and is adminis- 
tered by a faculty committee composed of repre- 
sentatives from these departments under the chair- 
manship of Dr. R. W. Moulton, Executive Officer, 
Department of Chemical Engineering. The follow- 
ing courses are now offered. 


Introduction to Nuclear Engineering 

Nuclear Instruments 

Nuclear Reactor Theory 

Nuclear Reactor Theory (Advanced) 

Nuclear Engineering 

Heat Transfer 

Nuclear Power Plants 

Tracer Techniques in Mechanical Engineering 
Measurements 

Nuclear Metallurgy 

Nuclear Metallurgy Laboratory 

Control of Radioactive Wastes 

Radiochemistry 

Radiochemistry Laboratory 
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THE APPLICATION OF DIGITAL COMPUTERS IN 
POWER-SYSTEM HYDRO-STORAGE PROBLEMS 


Davip L. JOHNSON 
Assistant Professor of Electrical Engineering 


Determination of the most 
economical use of hydro 
storage, either for an indi- 
vidual power system or in 
power-pool operation, has 
long been a problem solved 
by an impressive array 
of trial-and-error solutions. 
Each trial solution has re- 
quired extensive hand cal- 
culation and use of numer- 
ous lengthy tables. <A 
considerable practical back- 
ground involving familiarity 
with the individual system is necessary in order that 
the number of separate solutions not become pro- 
hibitive. The usual method of solution of hydro- 
storage problems has been to start with an “educated 
guess” or engineering assumption; then by lengthy 
trial corrections, a solution is obtained that is hoped 
to be the optimum one under the prescribed con- 
ditions. 

Generally, this type of solution leaves much to be 
desired. It imposes considerable responsibility upon 
the engineer who is establishing the starting values, 
it requires extensive and repetitive calculation, and 
often gives little assurance that the final answer will 
be the most satisfactory one. 

In the spring of 1955, a cooperative study was 
undertaken by Seattle City Light and the University 
of Washington Department of Electrical Engineer- 
ing, one phase of which was concerned with the in- 
vestigation of the Seattle City Light hydro-usage 
problems by digital computer means. The computer 
facilities of both organizations at that time were IBM 
604 electronic calculating devices. 

The basic problem investigated has been a limited 
one; however, the methods and techniques devel- 
oped can well be applied to more complex systems. 


Problem Definition 

The particular problem involved is one concerning 
the hydro-storage of the combined Seattle and 
Tacoma systems in operation as a part of the North- 
west Power Pool under a Bonneville Power Adminis- 
tration (BPA) contract. Each utility is required to 
submit annually a “computed demand table” which 
considers estimated peak and average loads and gen- 
erating capabilities to determine monthly peak and 
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average power requirements from BPA for the forth- 
coming year. The calculations to fix the generating 
capabilities must assume minimum water-year condi- 
tions. The problem thus becomes one of accurate 
estimated loads and economical hydro usage. This 
discussion will not be concerned with load estimates, 
although this is an area in which digital computers 
could easily be used to advantage. 

The contract is so set up that a utility may not be 
credited for covering a pessimistic load estimate from 
its own facilities even though it may have the hydro 
capability to do so. The contract cost of BPA power 
to each utility will be dependent upon the maximum 
power demand of that utility, either peak or average, 
as shown by the computed demand table, which then 
becomes a theoretical device based upon actual esti- 
mated loads and a standard minimum water year. 
The accuracy of such a study and the requirement 
that the system use the minimum water-year capa- 
bilities to the optimum is a problem of extreme eco- 
nomic significance during a given contract year. 
Penalty clauses in the contract consider the possibility 
of underestimating the maximum power required of 
BPA by the utility. 

At present in the combined Seattle-Tacoma sys- 
tems, there is only one major storage reservoir which 
is used economically for variable-head storage. The 
other reservoirs are primarily suited for maintaining 
nearly a fixed head for maximum peaking capability. 
The variable Ross reservoir, however, feeds not only 
into the facilities at Ross plant, but downstream on 
the Skagit River to two other generating facilities, 


OTHER BPA 
| RESOURCES 
L 
LOAD 
we 4 
/ 
; 
ROSS 
/ 
DIABLO 
/ 
GORGE 


Fic. 1. SCHEMATIC REPRESENTATION OF THE SKAGIT 
Hypro SYSTEM 
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Diablo and Gorge plants. A schematic representation 
of the Skagit hydro system is shown in Fig. 1. 
Over a monthly period the heads at the down- 
stream reservoirs remain essentially constant. The 
problem breaks down into one of a single variable: 
storage reservoir. It should be noted that the gener- 
ating facilities on the Skagit River supply roughly 
one-half of the total combined system generation. 
As the cost per megawait-hour of BPA power is 
dependent upon the use factor and therefore is a 
function of the maximum power deficit of the util- 
ities as shown by the computed demand table, it is 
most economical to make the maximum peak power 
deficit equal to the maximum average energy deficit 
during the month in which this deficit occurs. Keep- 
ing the monthly average energy deficit equal to the 
maximum value during the withdrawal period in the 
months preceding the month of the maximum peak 
deficit will maintain maximum storage and therefore 
permit the lowest possible value of maximum peak 
and average power deficit during the contract year. 


Means of Solution 

The methods of solution of the problem previously 
used by Seattle City Light engineers was one set 
by the BPA. The general criteria previously men- 
tioned were followed in a repetitive trial-and-error 
solution. In light of experience, certain withdrawals 
were used in the early months of contract calculation. 
(The contract period extends from July through 
the following July.) An estimated value of the 
maximum deficit and the month in which this deficit 
would occur was assumed. Calculation was then 
made to determine the storage required from Ross 
reservoir to meet this maximum value. When the 
calculations were complete to the assumed month of 
the maximum deficit, and the average energy deficit 
and peak power deficit were neither equal to each 
other nor equal to the assumed value, a correction 
was made upon the assumed value and the sequence 
was recalculated. Each of the numerous calculations 
required the use of extensive tables and six-place 
accuracy. After the correct results were finally ob- 
tained to the month of the maximum, which often 
falls in January, the process was generally repeated 
for successive months under certain limits until the 
final result was complete refill of Ross reservoir by 
the end of the contract period. Such a method re- 
quired numerous man-hours and was seldom free 
from the consideration that a numerical error may 
have invalidated the whole calculation. 

When this problem was approached from the 
standpoint of a digital computer solution, two gen- 
eral courses were possible. One was simply to 
mechanize the required calculation by use of the 
computer; the other was to investigate the basic 
theory of the over-all problem and try to improve the 
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method as well as the means of solution. The latter 
course was followed, and a simple method leading 
to a more exact solution was obtained. 

The value of the average power output of the 
Skagit system can be expressed as a function of 
Ross reservoir storage. The constants involved are 
the initial storage at the beginning of the period of 
calculation, the minimum water-year inflows of the 
Skagit averaged over the period of calculation, and 
the information relating the quantity of water re- 
leased at a given head to the quantity of power pro- 
duced. This last information had previously been 
taken from extensive tables relating head to H/K 
(megawatts per cubic foot-second). For use on the 
computer the method of least squares was used to 
obtain H/K as a function of reservoir storage. The 
results obtained by this method were very satisfactory 
—probably as accurate as the original data. This 
operation was required only at Ross plant, since the 
head and H/K were assumed to be fixed at the down- 
stream plants. 

The power output of the Skagit system as aver- 
aged over d days can be expressed as 


=Fp (H/K)p+Fe¢ (H/K)¢ 


[((H/K)g+(H/K)p]+M (1) 
, 


where F is the total plant flow and NF is the natural 
flow under minimum water-year conditions, Qy and 
Q, are Ross reservoir storage in second-foot-days, at 
the beginning and end of the period of calculation 
(d days), respectively, and H/K is the average H/K 
over the period, represented as a linear function of 
average reservoir storage. 

Assuming that the Skagit power is the only system 
variable, and that the capabilities of the rest of the 
system are fixed and known during each period, the 
average energy deficit for a given period can be 
determined : 


Prov Pav Pow tx 


= , (2) 
where L,, is the total average predicted load for the 
system over the time period, and P,, is the total 
average energy capability of the system over the 
time period broken down into the variable, Pay sx, 


* Subscripts Sk, R, D, and G refer to Skagit, Ross, Diablo, 
ani Gorge, respectively. 


and the fixed power capability of the system other 
than the Skagit power, P,y gx. 

The peaking capability of the system is assumed 
to be fixed and known for all of the facilities other 
than Ross plant. This is a valid assumption for 
normal operating conditions. The information re- 
garding peaking capability has previously been used 
in tabulated form relating peaking capability to 
head. As in the case of the Ross H/K factor, the 
method of least squares was used to determine the 
peaking capability of Ross plant as a function of 
reservoir storage. This was done in the form of a 
quadratic, 

R=AQhin BO min tC ’ (3) 
where Qin is the minimum storage in the Ross 
reservoir during the time period. 

The peak power deficit can now be expressed in 
the form: 

Peak det = Lpeak — Ppeak 

~ fix * peak R (4) 

—AQhin— 
where Lyeax is the total peak predicted load for the 
system over the time period, P,.4, is the total sys- 
tem peak capability over the time period, made up 
of the variable P,.., x and the fixed peaking capa- 
bility of the system other than Ross, Ppeax fx- 


Equations 2 and 4 are oversimplified in that they 
do not consider the average and peak power reserves 
allowed under contract agreement. In the case of 
the average power, the reserve is a function of the 
steam capability of the system and therefore can be 
incorporated into Pay six. The peak power reserves 
allowed are more complex in that they may be fixed 
by either a percentage of the total peaking capability 
or by the loss of the largest machine, whichever is 
the lesser. This problem involves double solutions 
of equations, in many instances using two sets of 
coefficients and a decision as to which should be 
accepted, but does not alter the general structure of 
the solution. In summary, it can be stated that both 
the peak power and average energy deficits can be 
expressed in terms of the Ross reservoir storage, 
average predicted loads, average and peak generat- 
ing capability of the combined systems other than 
the Skagit, and the average stream flows as fixed by 
the minimum water-year. This is done by means of 
quadratic equations which have fixed coefficients for 
any given period of time and given system. The 
value of Ross reservoir storage is the only variable. 


One of the criteria previously mentioned was that 
of maintaining the average energy deficit equal to the 
value of the maximum peak power deficit from the 
time of withdrawal to the period when the maximum 
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peak deficit is reached. This can be mathematically 
expressed by operating upon Eqs. 2 and 4 in the 
following way. The monthly peak power deficit, Eq. 
4, for the final month of the period under calculation 
is equated to the mean of the average energy deficits 
over the period starting with the beginning of the 
drawdown and terminating upon the final month. 
This means equating Eqs. 4 and 2, using coefficients 
of the final month in (4) and coefficients averaged 
over the period of calculation in (2). This operation 
will result in a single quadratic equation for each 
period starting with the month in which withdrawal 
is started and ending upon successive months until 
a maximum deficit value is determined. The result 
of this operation will be the minimum peak power 
deficit for the annual period and the value of the 
average energy deficit that is required from with- 
drawal to the month of the maximum peak deficit to 
hold this value at its minimum. 

When the average energy deficit to be used over 
this period has been obtained, monthly calculations 
are made, using Eq. 2 to determine the Ross storage 
required at the end of each month through the month 
of the maximum peak deficit, and Eq. 4 to give the 
resultant peak power deficit for each month. A simple 
routine to determine whether Qo or Q, is equal to 
Omin is used before Eq. 4 can be run. 

During the periods immediately following the 
month of peak deficit, the loads are still high and the 
inflows are low, so the peak and average deficits are 
again equated on a monthly instead of an average 
basis, and the storage required to give these deficit 
values is calculated by using a quadratic equation 
made up by combining Eqs. 2 and 4. This storage is 
then used in the peaking Eq. 4 for determination of 
the value of the peak and average deficits. 

After this period, the primary consideration is that 
of reaching total refill of the Ross reservoir before 
the beginning of the next annual withdrawal period. 
The determination of which method should be used 
to accomplish the refill is subject to many limitations. 
The reservoir cannot be withdrawn beyond a certain 
critical level, the deficits cannot at any time exceed 
the maximum, and a certain minimal flow must be 
maintained at all times. These criteria require many 
decisions and the possibility of several recalculations. 
For instance, if the critical level of the reservoir is 
exceeded or the minimal stream flow is not main- 
tained under one method of withdrawal or refill fol- 
lowing the month of the maximum, another modified 
method must be used. These choices can be, and of 
course for computer application must be, system- 
atized in such a way that one negating result suggests 
another means of solution to the computer. 

A sample tabulation of the per-unit results ob- 
tained by digital solution of the computed demand 
problem is shown in Table I. 


THE TREND IN ENGINEERING 


TABLE I 


PER-UN:T RESULTS OBTAINED BY DIGITAL SOLUTION 
OF THE COMPUTED DEMAND PROBLEM 


Skagit Combined 
Month-End Capability System 
Month Ross Power Deficit 
Reservoir 
Storage Peak | Average| Peak | Average 
July 1.000 2.300 |} 1.009 |-0.555 | 0.318 
August 1.000 2.300 | 0.633 |-0.374 | 0.905 


September | 0.953 2.271 | 0.709 |-0.141 | 1.000 


October 0.888 | 2.224] 0.751 | 0.287 | 1.000 
November | 0.780 | 2.141 | 0.921 | 0.681 | 1.000 
December | 0.695 | 2.078] 0.877 | 0.791 | 1.000 
January 0.556 | 1.956] 1.059; 1.000] 1.000 
February 0.421 | 1.818 | 1.011 | 0.988] 0.988 
March 0.361 | 1.744 | 0.732 | 0.889] 0.899 
April 0.310 | 1.653} 0.972 | 0.645 | 0.437 
May 0.493 | 1.653} 0.921 | 0.385 | 0.437 
June 0.924 | 1.897] 0.858 | 0.173 | 0.437 
July 1.000 | 2.244] 1.161 |-0.267| 0.437 


Computer Application 

The method of solution outlined above was in- 
itially programmed for the type 604 IBM computer. 
The complete solution of the sequence over the con- 
tract year was accomplished by the solution of from 
15 to 20 quadratic equations and 13 to 15 equations 
involving only simple substitution. It should be 
recognized that coefficients must be precalculated for 
each of the monthly peaking Eq. 4’s; for the monthly 
quadratics from which storage is calculated, Eq. 2; 
and for the quadratics which result from equating 
Eqs. 2 and 4. These coefficients, however, once cal- 
culated, will remain constant as long as the storage 
facilities of the system are unchanged. 

Calculation of the problem on the IBM 604 was 
accomplished in two sections, using a single set of 
computer panels. The two sections are specified as 
the “data-processing” section and the “operation” 
section. The data-processing section is a single run 
in which the system data are put into a form suitable 
for solution of the various equations and punched 
into the operation cards, which have previously been 
partially punched with precalculated fixed coefficients. 
In order to have consistent monthly programs, an 
operation card is prepared for each possible type of 
solution for every month. This means that, in most 
months, only two or three out of six operation cards 
will actually be used in the operation section. How- 
ever, it is not known before the operation run is 
started which of the various operation cards will 
be used for a given month. It should be emphasized 
that the only variable information to be punched on 
the data-processing cards will be that of estimated 
loads and outside system capacity. None of this need 
be precalculated. 
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Following the data-processing run, an operation 
run is made to determine the month and value of the 
maximum deficit. The maximum deficit value is then 
gang-punched into the monthly operation cards 
through the month of the maximum, and a run is 
made to fix the monthly storage used and the peak 
power deficit through this period. If the month of the 
maximum deficit falls in January, the February and 
March run is made by equating the average energy 
deficit and the peak power deficit in order to deter- 
mine this value and the storage used for these 
months. A calculation is next made to determine 
the constant average energy deficit required to give 
refill during the April-to-July period. This calcula- 
tion is in the form of Eq. 4. By using this value of 
average deficit, the monthly values of storage used 
(negative in this case) and peak power are deter- 
mined for each month. 


Unusual load schedules can create several vari- 
ances in this solution. If the method summarized 
should require more storage withdrawal than is 
available during April, a new average deficit must be 
calculated for that month in order to stay within 
system limitations. This equation again has the form 
of Eq. 4. This value might possibly be greater than 
the maximum deficit value, in which case the solution 
must be modified from February through April. This 
may be done by finding a constant average energy 
deficit for this period to result in complete withdrawal 
by the end of April. It is also possible that unusual 
load conditions might occur which would result in 
the February peak and average deficits being rela- 
tively low, with the March peak and average deficits 
higher than the maximum. In instances such as this, 
a constant average energy deficit which is below the 
maximum is determined and held through these two 
months. Provision is made for such possibilities in 
the program, and operation cards are punched and 
processed in the initial data-processing run to be 
prepared for these eventualities. 


The total solution of the problem for the contract 
year using the IBM 604 requires slightly less than 
one hour. The operating instructions are such that 
no familiarity with the problem is required by the 
operator. The average error for the year tested was 
about 0.05 per cent, with a maximum error of less 
than 1.50 per cent. When a number of cases of the 
type discussed had been run, it became increasingly 
evident that, while the 604 offered a practical solution 
to the problem, the use of this machine required 
extensive card handling, value comparisons and 
evaluations by the operator, and several operator 
decisions. These decisions were not those of a tech- 
nical nature, but of the order of “If A is bigger than 
B, do this; if B is bigger than A, do this; if A 
equals B, do this.” The 604 is certainly not the most 
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efficient machine by which to obtain a solution to a 
problem of this type. Therefore a solution using the 
larger and faster IBM 650 was considered. 


The IBM 650 with its 2000 words of storage and 
sizable order code is of such a magnitude that the 
hydro-storage problem can be solved from start to 
finish, using a stored program with no operator 
intervention. Because of the fact that the computer 
must make its own decisions and comparisons and 
chart its own course through the entire problem 
solution, numerous variations in the application of 
the general solution method are required. These 
variations, however, do not affect the general solu- 
tion method as it has been discussed, but only the 
means of presenting this method to the computer. 


The computer program is so set up that the system 
variables are supplied to the machine as a part of 
the input, and the machine is started through the 
program. The device itself will determine the course 
of calculations to be taken from each preceding oper- 
ation, evaluate the results, and stop only when an 
acceptable annual computed demand table has been 
obtained and punched out from the computer. At 
the time of this writing, the entire program has been 
organized and the section comparable to the opera- 
tion run on the 604 has been smoothed out and is 
in operation. From the results of this section of the 
over-all program it can be predicted that the total 
elapsed time required for a complete solution of the 
annual hydro-storage problem of the system studied 
will be approximately ten minutes. The accuracy will 
be slightly better than that available by means of 
the IBM 604. The limitation upon the accuracy is 
not one of possible significance in calculation, but 
that imposed by the mathematical representation of 
the tabulated material. If there were any necessity 
for higher ranges, the functional mathematical repre- 
sentations could be carried to greater accuracy or the 
table-lookup potential of the 650 could be used. 
These possibilities will be investigated if it appears 
that an improvement in accuracy of over-all results 
would be advantageous. Since much of the input 
information is only significant to three places, it 
seems doubtful that there is any merit in obtaining 
results to more than five places of accuracy. 


It must be recognized that the particular problem 
under study is a limited one. The computed demand 
table does not specify actual system drawdown pro- 
cedure, but is made under contractual specifications 
for the express purpose of satisfying certain contract 
requirements. The method used, however, could 
easily be extended to include actual operating condi- 
ditions and a greater number of variable and related 
reservoirs. Familiar computer methods for the solu- 
tion of simultaneous equations might well be incor- 
porated into the more complex solution. Addition of 
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certain specifications and limitations could also make 
the program adaptable for actual drawdown solution. 
In instances where tailwater elevation is a significant 
factor, this problem could be incorporated into Eqs. 
2 and 4 without changing the basic structure of the 
solution. 


As the problem is extended in scope, the solution 
will become considerably more complicated. The 
IBM 650 solution time is divided between actual 
calculation and decision operations; of these two 
sections, the decisions and the attendant changes in 
program take up the greater time. With the extended 
problem, the calculations will increase, for the solu- 
tion of simultaneous equations will be required; but 
the number of decisions will rise very sharply with 
the resultant increase in time requirements. It can 
be estimated that the time required by a computer 
solution of a hydro-storage problem may rise with 
the number of variable reservoirs between a linear 
and square relationship. 


The means of rapid calculation of withdrawal prob- 
lems gives rise to a number of studies concerned with 
things other than contract calculation. Mass studies 
of the effect of weather conditions might be of con- 
siderable service in evaluation of the effect of alter- 
nate load estimates such as those used in the contract 
calculations. Many changes in system require only 
minor changes in the precalculated coefficients, and 
the proposed method could be an integral part in 
analysis of future planning. Extensions of the method 
could be used to obtain exact solutions for control 
of large hydro networks or combinations thereof. 


Conclusion 


The solution of the problem involving economical 
use of hydro storage is one that lends itself well to 
digital computer solution. The preparation of any 
solution for computer application should be more 
than the mere mechanization of methods previously 
used, but must require a re-examination of the basic 
problem in terms of the most elementary mathematics 
and theory. This fact is well borne out in the ease 
of solution as demonstrated in this particular prob- 
lem. The* final method was certainly not the one 
with which the investigator initially approached the 
problem. As familiarity with the problem developed, 
however, simplifications became apparent, and the 
final method of exact solution gradually evolved. 
Because of the nature of the problem and the tre- 
mendous number of solution limits involved, the 
application of the general method of solution to a 
specific computer solution requires a frustratingly 
great number of computer decisions and variations. 
This can be handled moderately well, however, on a 
medium-sized digital computer such as the IBM 650 
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Approximate Laminar Skin Friction in Compressible Flows 
with Favorable Pressure Gradients* 


Ety S. LEVINSKY 
Graduate Student in Aeronautical Engineeringt 


Introduction 
Laminar skin friction 
with zero pressure gradients 
has been the subject of 
many papers. The classical 
solution was obtained by 
Pohlhausen! for incompres- 
sible flow by representing 
the velocity profile in the 
boundary layer by a power 
series in y/é, in which a 
E. S. Levinsky parameter is used to ac- 
count for the effects of 
pressure gradient. The accuracy of Pohlhausen’s 
solutions was shown by Thwaites? to be surprisingly 
good for favorable pressure gradients; however, for 
adverse pressure gradients, where the possibility of 
separation exists, the agreement was poor. Thwaites 
modified Pohlhausen’s procedure, improving its 
accuracy for the latter case. However, even this 
improved method failed to accurately predict 
separation. The present procedure is based on 
Thwaites’ method, and therefore is valid only for 

favorable pressure gradients. 

The calculation of compressible skin friction with 
pressure gradient was greatly simplified when, in 


-1949, Illingworth? and Stewartson‘* introduced a 


transformation of coordinates which reduced the 
compressible flow equations to the incompressible 
form, thus enabling incompressible methods to be 
used in obtaining a solution. Rott and Crabtree® in 
1952 used this transformation, along with the 
method of Thwaites, to calculate compressible skin 
friction. 

In 1953, Rott® generalized the Illingworth- 
Stewartson transformation, making it valid for any 
constant Prandtl number near unity. Following a 
procedure similar to that of Rott and Crabtree, this 
generalized transformation is used in the present 

* Arranged by the author from his M.S. thesis of the same 
title and his paper which won the regional award given by 
the Institute of the Aeronautical Sciences, West Coast Con- 
ference, Graduate Division, 1956, published by the IAS. 

tThe author was employed by Sete Airplane Company 
at the time he obtained his degree of M.S. in Aeronautical 


Engineering, 1956. He is now a senior aerodynamics engineer 
with the San Diego Division of Convair. 
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method of calculating skin friction. The viscosity 
law, wo/u=C(T>/T), first used by Chapman and 
Rubesin,’ is employed, rather than the more approxi- 
mate relation, po/u=7)/T, used by Rott and 
Crabtree. 
Many of the derivations presented in the discus- 
sion are in a highly abbreviated form.’ 
Nomenclature 
a=velocity of sound 
_ dU/U, 
a’=velocity profile parameter=— V/A 
B,=parameter defined in Eq. 8 
c=airfoil chord length 
C=constant in the linear viscosity law defined 
by Eq. 7b 
¢j,=local skin friction coefficient based on local 
free stream dynamic pressure=7/} 
cy~=local skin friction coefficient based on undis- 
turbed free stream dynamic pressure = 7/3 p,,u?, 
Cr..=average skin friction coefficient based on un- 
disturbed free stream dynamic pressure, de- 
fined by Eq. 27 
Cp=specific heat at constant pressure 
G=function of a’, defined by Eq. 17 
H=ratio of transformed displacement to momen- 
tum thickness, defined by Eq. 9 
J =function of a’, defined by Eq. 19 
k=thermal conductivity 
K=function of a’, defined hy Eq. 9 
1=function of m, defined by Eq. 12 
L=function of m, defined by Eq. 14 
m=boundary layer parameter, defined by Eq. 13 
M=Mach number=u/a 
p=pressure 
Q=function of a’, defined by Eq. 18 
R=gas constant 
Re,= Reynolds number based on local free stream 
conditions and distance x 
Re,..= Reynolds number based on undisturbed free 
stream conditions and length c 
=temperature 
u=velocity component parallel to surface 
v=velocity component normal to surface 
U,V=velocity components in transformed flow, 
defined by Eq. 5b 
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x =curvilinear coordinate along surface 
y=coordinate normal to surface 
X,Y =coordinates in transformed flow, defined by 
Eq. 5a 
 =ratio of specific heats 
6=boundary layer thickness 
A=boundary layer thickness in transformed flow 
A*=displacement thickness in transformed flow, 
defined by Eq. 11 
n=small positive constant 
é=momentum thickness in p5ysical coordinates 
©=momentum thickness in transformed flow, 
defined by Eq. 10 
=function of U/U,, defined by Eq. 20 
\=Pohlhausen boundary layer parameter 
= (8#/r) (du;/dx) 
u=absolute viscosity 
v=kinematic viscosity = u/p 
—=function of Prandtl number=1—+/¢ 
p=density 
o= Prandtl number = uC p/k 
r=shear force per unit area 
#=function of velocity 
Q=momentum thickness parameter, defined by 
Eq. 22 
Subscripts 
2 =condition at outer edge of boundary layer 
= undisturbed free stream conditions 
onditions on surface 
=reference condition in free stream 
c=reference length based on chord 
FP = flat plate 
Discussion 
We shall begin with the Prandtl boundary layer 
equations for two-dimensional laminar flow. They 
are the equation of continuity, the momentum 
equations in the x and y directions, the energy 
equation, and the equation of state, which are listed 
in this order as 


(pu) (0) = 0, (1a) 

=0, (1c) 

(1e) 


The energy equation (1d) may be considerably sim- 
plified by taking the Prandtl number o=uCp/k=1. 
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By multiplying the momentum equation (1b) by u 
and adding the result to (1d), and then making the 
additional restriction of zero heat transfer, the 
energy equation reduces to the same form as in 
non-viscous flow, 


T 
Ty 


Tani® assumed that a relation somewhat similar to 
Eq. 2 might be used to replace the energy equation, 
even when the Prandtl number was not unity. He 
used the 


1 uj 


where ee and @ is a function of u/m only, 
which must have the values, (0) =#(1)=0, if the 
temperature is to satisfy the boundary conditions. 
Equation 3 will be used in the present analysis, 
restricting the results to flows with zero heat trans- 
fer, but with arbitrary pressure gradients. 

The function @ was obtained by Tani, who repre- 
sented it by a power series in u/u, and then solved 
for the coefficients. He also stated the following 
integral relation which ® must satisfy, 


0 0 


Modified Illingworth-Stewartson Transformation 


For incompressible flow, the boundary layer equa- 
tions may be solved by integrating the momentum 
equation across the boundary layer. Such an in- 
tegral procedure would be difficult to apply to the 
compressible case with the boundary layer equations 
in their present form, because of the appearance of 
p, T, and yp, in addition to u and », as the basic 
flow variables. 

Illingworth? and Stewartson‘ introduced a trans- 
formation which reduced the continuity and momen- 
tum equations to the same form as for incompres- 
sible flow. Their transformation was limited to a 
Prandtl number of one. Rott’s® modified Illingworth- 
Stewartson transformation, which is valid for any 
Prandtl number, will be used in this paper, increas- 
ing the generality of the results. His transforma- 
tion may be written 


Po 


ox To OX’ Oy To 
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where x and y are the coordinates in the original 
flow, and X and Y are the corresponding coordinates 
in the transformed flow. The subscript 0 refers to a 
reference point in the free stream, and therefore the 


1 
isentropic relation, p:/p.=(71/To)7—', holds outside 
the boundary layer. 
The continuity equation (1a) may be expressed in 
a form similar to that for incompressible flow by 
applying the transformation equations, (5a) and (5b). 
Thus, in the transformed coordinates it becomes 


, aVv 
ox toy 7°: (6) 


The momentum equation in the x direction (1b) may 
also be transformed in a similar manner. First, 
however, it is necessary to introduce the viscosity- 
temperature law, 

Mo To 

7 

>: (7a) 
which was used by Chapman and Rubesin’ in their 
work on flat plates, where the variable C is defined 
by the highly accurate Sutherland law, 


T+216 


Now, by utilizing Eqs. 5a, 5b, and 7a, the momentum 
equation becomes, after some manipulation, 


U au Th 2 a? U; aU, 


Substituting for 7/7; from the revised energy 
equation (3), this reduces to 


aU dU, C (8) 


1 


Equation 8 is the momentum equation in the trans- 
formed flow. It is of the same form as given by 
Tani.’ The variable temperature, density, and vis- 
cosity which appeared in the original form (1b) have 
been eliminated. Except for the appearance of @ and 
the deviation of C from unity, Eq. 8 is identical with 
the incompressible momentum equation. 


where B,=(y—1) (M?/2) 
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Integration of the Momentum Equation 


The solution of Eqs. 6 and 8 can be accomplished 
by formally integrating Eq. 8 across the boundary 
layer, the relationship expressed by the continuity 
equation (6) being used in the process. The resulting 
integrals are then evaluated by representing the 
velocity U by a power series in Y/A, where A is the 
boundary layer thickness in the transformed flow. 
Because of the appearance of @ and C in the momen- 
tum equation, formal integrals involving these func- 
tions will also occur. Tani evaluated these integrals 
by representing @ and C also by power series in 
Y/A, and some of his results will be employed in 
the following analysis. 


By integrating Eq. 8 from Y=0 to Y=A, setting 
the reference temperature 7)>=T,, and making use 
of the continuity equation (6), the momentum 
equation becomes 


Vo aX 


which is identical to the integrated momentum 
equation in incompressible flow, except for the term 


in §B,K, where K is given by Tani as K= =t foay 
, is: de- 


= ay, (10) 


H is given by the ratio, H=A*/@, where A*, the 
transformed displacement thickness, is 


ay, | (11) 


while / is obtained di the equation, 


wf (12) 
U; oY 


Now, for incompressible flow, Thwaites? expressed 
both H and | in the terms of the parameter m, where 
dU, 
Vo “ax 
The determination of H and / in terms of m will not 
be covered in detail, except to remark that both 
functions were evaluated by Thwaites independently 
of Eq. 9, and therefore his results are valid for the 
present analysis, even though Eq. 9 is not identical 
to the incompressible equation. By writing Eq. 9 
in the form, 


ax j (14) 


The transformed momentum thickness, @ 
fined as 


m=— 


(13) 
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it is seen that, for incompressible flow, with B,=0, 
the right hand side, L(m), is a function of m only. 
Both H(m) and l(m) have been plotted in Fig. 1, so 
that L(m) may be obtained from Eq. 14. Thwaites 
approximated L(m) by the linear function, L(m) 


=0.45+6m, and thus succeeded in integrating 
Eq. 14 directly to 
2 xX 
@”=0.45 f UsdX. (15) 
Vo 0 


Equation 15 may be applied to compressible flow 
only for a Prandtl number equal to one, for which 
£=:(. Some estimate of the error involved in using 
this equation for a Prandtl number of 0.72, corre- 
sponding to air at standard conditions, and for 
intermediate supersonic Mach numbers can be found 
from the consideration that, for these conditions, 
—£ B,K is approximately 10 per cent of 2+H(m). 
Such an error will be considered too large to ignore 
in the development of Eq. 15 from Eq. 14, and 
therefore Eq. 15 will not be used in the subsequent 
analysis. 

A direct relation between @ and X without the 
term in K may be obtained by integrating the 
product of the transformed momentum equation (8) 
with U across the boundary layer. By making use 
of the transformed continuity equation (6) and the 
relation given by Eq. 4, the result is 


Or +L (14-3658) St 
16 
where 
=e*/o, (17) 
Yer, (18) 


ands J= (1-1) ( (19) 


and where the quantity ©* is defined as 


4U U? 
@*= U, (20) 


Tani evaluated G, Q, and / for both adverse and 
favorable pressure gradients by representing U’ and 
C as a power series in Y/A, and obtaining the 
constants in the series from the boundary con- 
ditions at Y/A=0 and Y/A=1. Tani chose to 
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solve for the constants in terms of the function, 
a’(X)=0(U/U,)/a(Y/A), rather than for m as is 
usually done. 


Momentum Thickness 


The momentum thickness in the transformed flow 
may be obtained by integrating Eq. 16 with respect 
to X. The additional assumptions were made that 
the variables, Q/G? and /, could be replaced with 
the constants, 0.44 and 0.184, respectively, corre- 
sponding to their values for zero pressure gradient; 
and also that the variables C; and B, could be ap- 
proximated by their free stream values, C, and B,. 
Equation 16 then integrates to 


2 
1+0.184( -1)] x 
Vo 


x 
gx on) 


which is remarkably similar to Eq. 15 obtained by 
Thwaites. Although Eq. 21 in its present form could 
be used to calculate skin friction, it is more con- 
venient to deal with quantities in the original com- 
pressible flow, rather than in the transformed flow. 
If we define the momentum thickness @ in the 
original flow as 0= dy, and ex- 
Pit; 

press the other seinatinies in Eq. 21 in terms of the 
original flow variables, Eq. 21 becomes 
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(22) 
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where the local Mach number, M;, rather than 
velocity ™“, has been used for convenience. The 
quantity Q is referred to as the momentum thickness 
parameter. Its evaluation over the entire surface is 
a prerequisite for determining skin friction. 
Skin Friction 

Skin friction on a surface results from shear forces 
which are proportional to the velocity gradient in 
the boundary layer normal to the surface, the con- 
stant of proportionality being the absolute vis- 
cosity, uw. Thus, the shear force per unit area, 7, at 
any point along the surface may be expressed as 
raul | . By making use of the definition of / 
(Eq. 12) and by writing 7 in terms of variables in 
the original flow only, this becomes 


Ti 
( T. ). (23) 


where the reference temperature, 7»), and hence 
viscosity, uo, are equal to the values at the surface. 
Then, by defining the local skin friction coefficient, 
cy,, based on the local dynamic pressure, p:u?/2, as 


2 
- , Eq. 23 becomes 


Q M, 


w 


where the local Reynolds number Re,= Rees, 
Q is given by Eq. 22, and / may be found from Fig. 1 
after m is known. 

We may find m in terms of quantities in the trans- 
formed flow from Eq. 13. However, the calculation 
of skin friction is greatly simplified if m is also 
written as a function of the original flow variables. 
Equation 13 then becomes. 

Ty y dM,/d(x/c) 
The local skin-friction coefficient on the surface may 
be obtained by evaluating Eqs. 22, 25, and 24 in 
that order, with / being found from Fig. 1. Before 
these equations may be evaluated, the Mach number 
M, and temperature 7, in the local free stream must 
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be obtained, the initial conditions being the free 
stream Mach number, M,, and temperature, T,,, 
and the surface or wall temperature, T,. 

Up to this point, only the local skin friction at 
each point on the surface has been considered. The 
average skin friction coefficient, Cr, based on the 
undisturbed free stream dynamic pressure, p,.u2,/2, 
is of practical importance, for it leads to the total 
viscous drag and may be readily compared with flat- 
plate values. By defining cy, as the local skin- 
friction coefficient based on free stream dynamic 
pressure, it can be readily shown that 


CAV Re, 


where Re... is the free stream Reynolds number 

based on the chord length, c, and is defined as 

Reo, = . The average skin-friction coefficient, 
Hoo 

Cr, is related to cs, by the expression, 


(27) 


which could be used, together with Eq. 26, to obtain 
the average skin-friction coefficients on two-dimen- 
sional shapes, if ‘t were not for the infinity which 
occurs in the integrand for x/c=0. This difficulty is 
overcome by assuming that the integrand varies as 
(x/c)* in the interval 0=x/cZn, being a small 
positive number such that 0<4<<1. (For the bicon- 
vex airfoil calculations, 7=0.01.) Equation 27 is 
then replaced by the equivalent relation, 


Equations 26 and 28 enable Cr,~/ Re, to be ob- 
tained, after c;,~/ Re; has been evaluated over the 
entire surface. 


Results for Zero Pressure Gradient 


For the particular case of zero pressure gradient, 
the skin friction and momentum thickness reduce to 
very simple-expressions, and a comparison may be 
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FiG. 2 FLAT-PLATE SKIN FRICTION 


made with the results of other investigators. For 
this ease, m=0 and /=0.22 and, therefore, 


C 
=1. 9 
Ce. 1.3264) (29) 


which is identical with Eq. 58a of Chapman 
and Rubesin.?’ Equation 29 is compared in Figs. 
2a and 2b, for 7,,=648° R and 72° R, respectively, 
with the method of Rott and Crabtree,® and with 
the calculations of Young and Kirkby.’ The highly 
accurate but involved calculations of Crocco" have 
also been included for comparison. Equation 29 
agrees reasonably well with Crocco’s results. Rott 
and Crabtree’s method, however, indicates a con- 
stant value of Cr,~/Ke., with Mach number and 
temperature. This results from their use of the 
viscosity law, uo/u=7./T, which is equivalent to 
the Chapman-Rubesin relation used in the present 
paper, with the constant C, =1 for all conditions of 
temperature and Mach number. Thus, by setting 
C,=1, they eliminated the Mach number- and 
temperature-dependent terms from their theory. 
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Fic. 4. LOCAL SKIN-FRICTION COEFFICIENTS WITH Co=1 
FOR FIVE-PERCENT PARABOLIC ARC AIRFOILS 


Young and Kirkby’s results, based on the viscosity 
law, po/u=(To/T)**, indicate that Cr,~W/ Re. is 
independent of T., and thus also fail to agree with 
Crocco’s calculations. On the basis of these com- 
parisons for zero pressure gradient, it was concluded 
that to correctly predict the variation of skin fric- 
tion with Mach number and temperature, a viscosity 
law of the type defined by Eqs. 7a and 7b must be 
employed, even when the pressure gradient is 
not zero. 
Airfoil Skin Friction 

Skin friction calculations, by the present method, 
were then made for a 5-percent-thick parabolic arc 
airfoil section. The detailed calculations are given 
in Reference 8, but some of the results have been 
included in Figs. 3 through 6. The computations 
were carried out for M,=1.5 and 4.0, T7,=136° R 
and 392° R, and «=0.72 and 1.0. 

The variation of with x/c for M, = 4.0 and 
o=1.0 is indicated in Fig. 3 for both temperatures. 
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Calculations by the method of Rott and Crabtree 
and the two points by Weil” are also shown for 
comparison. Weil took C,.=1, and it is therefore not 
surprising that his values agreed well with those by 
Rott and Crabtree’s method. 


By dividing the skin-friction coefficient by 
the results of the present method become nearly 
coincident with those of Weil and with those ob- 
tained by Rott and Crabtree’s method, as indicated 
in Fig. 4. Similar results were found for M,=1.5, 
and it therefore appears that, for favorable pressure 
gradients, the skin friction is nearly proportional to 
VC, , in agreement with the results of Chapman 
and Rubesin for zero pressure gradient (Eq. 29). 


Prandtl Number Effect 

The foregoing calculations were then repeated for 
a Prandtl number of 0.72, corresponding to that for 
air in this temperature range. The change in skin 
friction was found to be very small (Table I). It was 
therefore concluded that, at least for the Prandtl 
number range from 0.72 =o =1.0, this effect could be 
ignored. In addition, it is recommended that o=1 
be used in the calculations rather than the actual 
value, because this simplifies the computations 
considerably. 

TABLE I 


AVERAGE SKIN FRICTION SUMMARY 
Five-Percent Parabolic Arc Airfoil Section 


TEMPERATURE Che VRea, 
Ma o Flat {Increase 

Static Total | Airfoil| Plate | (%) 
13 1.636 1.28 28 
569° R | 1.636 1.28 28 
¢ 1.647 1.29 28 
13 569° R | 1.647 1.29 28 
4.0 {1.0 392°R 1.365 1.08 26 
569° R | 1.635 4.23 31 
AO 1.364 | 1.106 23 
569° R 1.633 | 1.266 29 


Average Skin Friction 

The ratio of the average airfoil skin-friction co- 
efficient, Cr,, to that for a flat plate, Crorp, at the 
same free stream Mach number and temperature, is 
plotted in Figs. 5a and 5b for o=0.72 and 1.0, re- 
spectively. The results of Young and Kirkby," 
which, as mentioned previously, are based on the 
temperature-viscosity relation yo/u=(To/T)*, have 
been included in Fig. 5a for comparison, while those 
of Boison," which are based on po/u=T7)/T, have 
been included in Fig. 5b for the same reason. Note 
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that, while the basic values of Cr,, for Boison’s and 
Young and Kirkby’s calculations differ from those 
of the present method, the ratio Cr,/Crarp is es- 
sentially the same. 

A summary of the results obtained by the present 
method is given in Table I. Here it is shown that, 
for the cases considered, the average skin friction for 
the 5-percent airfoil section is 23 per cent to 31 per 
cent greater than for a flat plate. 


Conclusion 


A method of estimating compressible, laminar, 
skin friction on two-dimensional bodies with favor- 
able pressure gradients has been presented. Equa- 
tions 22, 24, and 25 are required to obtain local skin- 
friction coefficients, while Eqs. 26 and 28 are needed 
for average skin-friction coefficients. 

(Continued on page 22) 
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NOUN PHRASE ANALYSIS USING THE IBM 650 COMPUTER 


ARISTOTELIS D. STATHACOPOULCS 
Graduate Student in Electrical Engineering 


At present, machine trans- 
lation is accomplished by a 
word-for-word translating 
process, including some id- 
iomatic phrases treated as 
independent units and trans- 
lated as such. It is easily 
apparent that translation 
of this kind carries the de- 
sired information, but it 
requires special effort by the 
reader in order to be under- 
stood. The output of such 
a translating device will 
appear in the source-language word order and it 
will include multiple translations for words with 
several meanings in the target language. In other 
words, such an output is not in conventional English 
and requires polishing and rewriting before it can 
compare with the output of a good human translator. 

Rearrangement of word order and selection of the 
contextual meaning of a word constitute two of the 
main problems of machine translation. A device 
which will overcome these difficulties should be able 
to perform certain logical operations which the mind 
of the human translator performs unconsciously. 

This fact suggests the use of a modern computer 
with permanent stored program to analyze the 
source-language text and perform the changes re- 
quired in the word-for-word equivalent target-lan- 
guage translation. 

One particular problem is recognition of the word 
series consisting of article, adjectives, and final noun, 
called a noun phrase, which replaces the simple noun 
in more complex sentence structure. The following 
discussion describes a method for recognizing and 
analyzing a noun phrase in any sentence. 

A noun phrase consists of articles, adjectives, and 
nouns. However, many of the constituents of a noun 
phrase may belong to more than one grammatical 
class. These words must be assigned to their correct 
class on the basis of context and translated accord- 
ingly. If this is to be done by a computer, each word 
must be accompanied by a class tag operated on by 
the computer program as it performs the desired 
logical operation while analyzing a group of words. 

Binary numbers were selected as tags and the 
intersection sum of logical algebra was selected as the 
main operation for the program. In taking the inter- 
section sum of two binary numbers the answer will 
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have a one whenever we add two ones in the same 
column, while it will have zeros in any other case. 
The example below illustrates what is meant by the 
intersection sum. 
11001 
<— binary numbers 


10010 


10000 <— intersection sum 
During the following discussion, by “‘addition” we 
mean the intersection sum. 
The tags selected for this analysis are the fol- 
lowing : 


110001 
adjective ..... .... 100001 
100000 


These tags have the following properties of logical 
algebra : 

art.(\art. =art. art.(\adj. =adj. art.(\noun = noun 

adj.(\adj.=adj. adj./\noun=noun 
noun/ \noun = noun 

The symbol () indicates the intersection sum. 

Each word will be represented by one tag for each 
grammatical class to which it belongs. Thus, words 
belonging to several grammatical classes will be rep- 
resented by a series of tags in a predetermined order. 
For example, the word “book” will have the follow- 
ing series of tags: 

100001 100000 111111 

The tag 111111 is the tag assigned to any verb. 

The following description outlines the essential 
features of the flow chart of the program as shown in 
Fig. 1, which provides for a noun phrase of up to six 
words. 

The first tag of the first word is examined to 
determine whether it marks the beginning of a noun 
phrase or not. If not, a further analysis follows to 
determing the nature of the tag. However, if the 
first tag indicates the beginning of a noun phrase, 
the routine continues as follows: The first tags of 
the next five words are added consecutively to this 
tag, and after every addition, the result is tested to 
determine whether we have a noun phrase or not. If 
the first tags fail to indicate a noun phrase, the second 
tags are examined until a noun phrase tag is formed. 

The following example of a three-word noun phrase 
will demonstrate the performance of the program: 

The Australian team appeared. 
110001 100001 100001 111111 
000000 000000 100000 000000 
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Initially, the first tag for the word “the” is an- 
alyzed, and since it is an article tag, it marks the 
beginning of a noun phrase. Following this, the first 
tag of the word “Australian” is added to the previous 
tag and the result is 100001 which is not the noun 
phrase tag. Consecutively, in the same fashion, the 
first tags for the words “team” and “appeared” are 
added, but still these tags fail to form a noun phrase 
tag, 100000. Then the first tag of the word “ap- 
peared” is replaced by the second tag and a new 
intersection operation is performed. Still no noun 
phrase is obtained. Then “appeared” is si: ste 
the process. Next the first tag of the word ‘ 
is replaced by its second tag and a new eels od 
sum is obtained. The final result is 100000, which 
indicates the end of the noun phrase. 

At this point, a card will be punched by the machine 
with the tags 110001, 100001, 100000, which are the 
correct tags of the analyzed noun phrase “The Austra- 
lian team.” In this process the noun phrase was 
identified and “team” was found to belong to the 
noun class in this context. 


Since the IBM 650 computer is a bi-quinary 
machine and is not capable of performing intersection 
summations, a special subroutine has been included 
in the program which is designed to transform the 
normal summation result to the corresponding inter- 
section sum. 

The program is able to determine the beginning 
and the end of the noun phrase and assign the cor- 
rect grammatical class to each constituent. It also 
demonstrates the possibility of applying mathematical 
operations to sentence analysis and the use of a com- 
puter to carry out the process. 

A complete program of this kind, constructed for 
analyzing sentences of the source language, could 
then be incorporated into a word-for-word translator 
and work on the tags of the source-language words. 
It could be expanded further to be capable of trans- 
forming the source-language analysis results to the 
corresponding target-language analysis and thus im- 
prove the accuracy, understandability, and readability 
of the final translation. 


Laminar Skin Friction in Compressible Flows 


(Continued trom page 19) 


Calculations carried out by this method for both 
airfoils and flat plates justified the following con- 
clusions: 


1. To accurately predict the variation of skin fric- 
tion with Mach number and temperature, a 
viscosity law of the type employed by Chapman 
and Rubesin (Eqs. 7a and 7b) must be used. 

2. Laminar skin friction is very nearly proportional 
to VYC,, for zero and favorable pressure 
gradients. 

3. The effect of Prandtl number on laminar skin 
friction is negligible. 

4. The ratio of airfoil to flat-plate skin friction is 
only slightly affected by the type of viscosity 
law used. 
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Zone Melting: A Physical Method for Controlling 
Impurities in Metals 


DouG tas H. POLONIs 
Assistant Professor of Metallurgical Engineering 


Significance of Impurities 

A thorough knowledge of 
the properties of pure metals 
is the ideal starting point 
for understanding the struc- 
ture and behavior of com- 
plex alloys. Trace impur- 
ities which are inevitably 
present in metals will mask 
the physical properties of 
the ultrapure material. Con- 
sequently, very little is 
known about the properties 
of ultrapure metals at the 
present time, owing to the difficulties in achieving 
high purity. 

The usual laboratory specimens which are specified 
as pure metals always contain trace impurities and 
are therefore not really pure in the true sense. Lead, 
for example, may be specified as 99.999-+ per cent 
pure but detectable amounts of bismuth, copper, iron, 
and other metals will also be present. Until recently 
the attainable purity in the common metals has been 
limited by the efficiency of chemical extraction and 
refining processes. For many metals the purity 
obtainable falls far short of that for lead; for exam- 
ple, titanium can be prepared only 99.9+ per cent 
pure by thermal decomposition of titanium tetra- 
iodide, and this is considered the highest grade titan- 
ium available. In order to increase the purity of 
metals far beyond that possible by chemical methods, 
the zone-melting technique has been developed and is 
becoming increasingly more important both in re- 
search and industry. 
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Nature of Impurities 

Residual impurities in metals are traceable to sev- 
eral sources : 

a) the ore from which the metal is obtained ; 

b) contaminants introduced during smelting, melt- 

ing, casting ; 

c) residuals from refining ; 

d) contamination during heat treatment. 

The increasing importance of the reactive metals 
in the aircraft industry and nuclear engineering has 
made item (d) more significant than in the past. 
The contamination of zirconium and titanium by 
small quantities of oxygen, nitrogen, carbon, and 
hydrogen has serious effects on physical properties. 
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Impurities may be present as dispersed inclusions 
such as oxides, sulfides, and silicates; as gases which 
cause internal cavities or blowholes; or as structural 
constitutents such as solid solutions or second phases, 
depending on the solubility of the impurities in the 
metal. Most minor trace elements are present in 
such small quantities that they form solid solutions 
with the metal; thus the impurity, in this case, is 
randomly distributed throughout the base metal or 
solvent on an atomic scale. These solid-solution im- 
purities are of particular concern in the present dis- 
cussion. 


Effects of Minor Impurities 

Trace impurities can have significant effects on the 
properties of metals. Foreign atoms in solid solution 
in a metal structure will cause some deviation of 
mechanical and electrical properties from those of 
the pure metal, particularly in view of the difference 
in atomic diameters and the probable valency differ- 
ence of the solute. It is interesting at this point to 
consider several familiar examples where trace 
elements exert a definite influence which is explain- 
able in terms of current theory of the metallic state. 


The effects of minor impurities on mechanical 
properties is well illustrated by the case of iron 
containing carbon concentrations ranging from 10-* 
to 10°? weight per cent. Tensile testing reveals an 
upper and lower yield-point phenomenon which is 
familiar to most engineers in stress-strain relations 
for low-carbon steels. This anomalous yielding effect 
has been attributed to the clustering of carbon atoms 
in the form of an “atmosphere” around dislocations, 
thereby increasing the stress to initiate dislocation 
movements to produce slip.'. Recent research has 
shown that some nonferrous metals, including molyb- 
denum and cadmium, also exhibit upper and lower 
yield points when they contain interstitial solutes. 


The development of semiconductors has presented 
some unusual metallurgical problems, particularly in 
connection with trace impurities in germanium and 
silicon and their effect on electrical properties. The 
distribution and types of small concentrations of for- 
eign atoms, including phosphorus, arsenic, and anti- 
mony, in these semimetals has a great influence on 
their usefulness in electronics, notably as transistor 
materials. As a result of efforts to produce extremely 
pure silicon and germanium, the zone-refining process 
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Fic. 1. Partial phase diagram of a binary system, A-B, 
showing limited solid solubility and lowering of liquidus 
temperature with increasing solute concentration. 
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FiG. 2. Schematic representation of interface temperature 
at successive stages of solidification and the corresponding 
solute concentration in solid and liquid at the interface. 
Note: Shaded areas represent excess solute accumulation in the liquid 


at successive interface positions during the early stages of solidifi- 
cation. 


was developed at the metallurgical laboratories of the 
Bell Telephone Company.** 

The zone-refining process is extremely simple in 
principle and, in this respect, is an elegant example 
of the application of fundamental principles in solv- 
ing practical problems. The process is now being 
applied to the purification of many metals, primarily 
for engineering research purposes. This paper will 
analyze the principles of controlled segregation and 
zone refining. 


Principles of Zone Melting 


Zone melting is similar in principle to the fractional 
crystallization process which has long been used as 
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a separation method by chemists. Briefly, the process 
involves purification by passing a molten zone along 
a bar of metal, thereby picking up impurities in the 
liquid and leaving behind a solid of lower impurity 
content. A basic requirement of this process is a 
higher solubility of impurities in the liquid phase than 
in the solid phase. 

Details of zone melting can be most easily under- 
stood by referring to a portion of a simple phase dia- 
gram of the type shown in Fig. 1. In the system 
shown, the solute B lowers the freezing temperature 
of pure component A and causes solidification over 
a range of temperatures in accordance with the Phase 
Rule for a two-component system containing two 
condensed phases. This is a hypothetical case since 
most metals will contain more than one impurity, and 
consequently a binary diagram is an extreme simpli- 
fication for actual conditions. It should be empha- 
sized that freezing of a liquid metal is a dynamic proc- 
ess; therefore, interpretation of freezing processes 
by reference to a phase diagram provides only an 
approximate picture of the actual conditions, even 
during low transformation rates. Nevertheless, it is 
possible to estimate conditions at a solid-liquid inter- 
face during solidification if one assumes that the 
solute concentrations at the interface for any temper- 
ature in the «+liquid region are approximately those 
predicted by the phase diagram. 


Consider the directional solidification of a long 
rod-like charge of molten alloy of composition Co 
in which the temperature is a minimum at the ex- 
treme left and increases uniformly along the charge. 
The first solid forms when the cool end of the charge 
becomes lower than temperature 7, and will have an 
approximate composition C, which is considerably 
lower than the average solute concentration of the 
original liquid. The formation of solid phase estab- 
lishes a solid-liquid interface which will move along 
the charge as solidification proceeds. This is shown 
schematically in Fig. 2, which illustrates the ideal 
case for uniform longitudinal heat flow. The solute 
depletion in the first solid to form provides the basis 
of purification. As freezing continues, the temper- 
ature at the interface must decrease and, accordingly, 
the composition of the solid phase forming at any 
temperature will be given approximately by the tem- 
perature intercept on the solidus line; i.e., the com- 
position follows: C,, C,, C,, Co between temperatures 
T, and T,. Eventually, composition of the solid at 
the interface reaches a maximum of Cy when the tem- 
perature reaches T7,. While the composition of the 
solid changes as described above, a simultaneous 
build-up of solute in the liquid takes place at the 
interface, where the composition follows the liquidus 
line between temperatures T, and 7,. It is seen that 
the maximum solute concentration in the liquid 
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reaches a limiting value C,; this concentration build- 
up in the liquid compensates for the depletion in the 
solid. If the molten charge is long enough so that a 
portion of the charge remains molten when the inter- 
face temperature reaches T,, then a condition of 
“kinetic equilibrium” is reached wherein the solute 
concentration at the interface in both solid and liquid 
(Co in solid and C, in liquid) and the interface tem- 
perature remain approximately constant as the inter- 
face moves along the charge. The last liquid to solid- 
ify will contain the excess solute which was removed 
from the first solid to form. Consequently, the im- 
purity concentration is built up in the final solid to 
freeze at the end of the rod. The concentration dis- 
tribution of impurity along the rod after solidification 
is completed is shown in Fig. 3. It is important to 
note the following facts regarding this process: 


1) Purification is effective only during the early 
stages of solidification, when the interface temper- 
ature is changing between 7, and 7,. During this 
stage the solid forming has a lower impurity content 
than the initial liquid, and some purification is 
achieved. This is called the transitional stage as op- 
posed to the steady-state stage, which is reached at 
temperature T, and continues until freezing is com- 
pleted. 

2) Diffusion is occurring in both the solid and 
liquid phases at all times. The amount of solute 
redistribution in the solid after freezing is negligible 
in practice. The transfer of solute across the inter- 
face and diffusion through the solid to level the con- 
centration gradient is so slow that it does not alter 
the solute distribution significantly during the time 
required to solidify the charge. There will be appre- 
ciable redistribution of solute in the liquid due to 
- diffusion and convection which result from the high 
atomic mobility in the liquid. 

3) The process described thus far is merely an 
idealized picture of solid-solution segregation com- 
monly referred to as “coring” and found in many 
commercial alloys in the as-cast condition. 

An important factor which indicates the efficiency 
of segregation in solid-solution freezing is the sep- 
aration of solidus and liquidus lines. This can be 
expressed as the Distribution Coefficient k at any 
temperature below the melting point for Fig. 1. Thus, 


Solute concentration in solid _— C Cs 
Solute concentration in liquid 3 


At temp T,, b=. At temp Te 

Co 
Reference to Fig. 1 reveals that the first solid to 
form has a concentration kCo. Consequently, in order 
to have a purification process, it is necessary that 
k<1. The purification effect ends when the solid 
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IMPURITY BUILD-UP 
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Fic. 3. Approximate concentration distribution along a rod 
after directional solidification. 


(a) 
STATIONARY HEATER 
(RESISTANCE OR INDUCTION) 


OF 
ROD 
MOLTEN 


MOVEMENT 
HEATERS 


Fic. 4. Schematic ram of zone-refining technique. Rod 
specimen is yg rough heater coil and molten zone 
moves along the specimen 

(a) Single pass setup. (b) Multiple zone-refining setup, 
with three molten zones passing simultaneously. 


reaches concentration Co, in which case the corre- 
sponding concentration at the interface in the liquid 
has a value Co/k. 


Single Pass Zone Refining 

The purification effect of the transitional solidifica- 
tion stage previously described can be successfully 
achieved by passing a molten zone along the speci- 
men. In this method only a fraction of the length 
of a specimen is molten at any time, so the impurity 
content builds up in the molten zone and is moved 
to the end of the specimen. The impurity content is 
thus concentrated in the last liquid to freeze. The 
concentration of impurity will build up in the liquid, 
with normal diffusion processes tending to level con- 
centration gradients in the liquid zone owing to rela- 
tively high atomic mobility in the liquid phase. The 
mobile molten zone can be most simply controlled by 
drawing a rod-like specimen through an induction 
coil or resistance furnace as shown in Fig. 4. In this 
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way only a localized portion of the charge will be 
heated above the melting point at any time. 

If a molten zone of length / is started at one end 
of a metal bar and passed slowly along its length, 
the segregation of impurity is similar to that previ- 
ously described for directional solidification. As the 
zone moves to the right, the first solid to freeze as 
the zone moves will have a low solute content, and 
the solute concentration will build up in the molten 
zone until the solute capacity of the liquid is reached. 
It is important to note the significance of the rate of 
movement of the zone. The ability of the liquid zone 
to contain impurity atoms depends largely on the 
diffusion of impurities through the liquid if the max- 
imum capacity of the liquid for impurity is to be 
utilized. Rapid zone movement should be avoided 
for efficient purification of the solid, since some time 
should be allowed for redistribution of the accumu- 
lated impurity within the molten zone. This does not 
present a serious experimental problem since normal 
convection effects in the liquid cause fairly rapid dis- 
tribution of impurity in the zone. After the zone has 
picked up its capacity of impurity, the solid freezing 
at any time during the remainder of the process will 
have an approximate composition Co. The excess 
solute deposited in the final region to solidify corre- 
sponds to the amount which was removed during the 
early stages and swept along in the molten zone. The 
solute distribution curve of Fig. 3 also applies approx- 
imately to the passage of a single molten zone along a 
specimen. 


Repeated Pass Zone Refining 

In order to achieve a high degree of purity along 
most of the length of a specimen, it is necessary to 
repeatedly pass molten zones along the specimen. 
Passage of successive molten zones progressively in- 
creases the purity of the initial regions at the starting 
end of the charge and extends the transitional, or 
purification, stage further along the bar with suc- 
cessive passes. Since with each successive pass there 
is less impurity to pick up in the initial end of the bar, 
the zone moves farther along the bar with each suc- 
cessive pass before reaching its capacity of impurity. 
Consequently, after a large number of passes in 
one direction, the impurity is concentrated at the 
final end of the specimen. The degree of impurity 
obtained throughout most of the bar is limited by the 
number of molten passes. The smaller the value of 
k=(Cg/C_,), the more efficient the process per pass. 

It should be noted that repeated pass zone melting 
can result in tapering a rod-shaped specimen. A 
slight taper may result even from a single pass and 
this effect will be quite noticeable after several passes. 
Producing Homogeneous Solute Distribution 

In some cases it has been desired to produce metal 
specimens containing dilute alloy concentrations in 
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SOLID CORNERS OR “FINS” 


\ 2-REVERSE PASS OF MOLTEN ZONE 


SOLUTE CONCENTRATION 


1-SINGLE PASS OF A MOLTEN ZONE 


2 
Ci 


ELEVATION 


DISTANCE ALONG SPECIMEN 
(Left): Fic. 5. Solute distribution along a specimen after 
single pass and reverse pass of molten zone. 


(Right): Fic. 6. “Cage” zone-refining technique which eli- 
minates ceramic containers. 


Note: Bar may be of polygonal cross section other than rectangular. 


which the solute is uniformly distributed throughout. 
This is difficult with normal casting because of solid- 
solution segregation. Uniform solute distribution can 
be achieved quite simply by the following three 
methods : 

a) Directional solidification: This method will pro- 
duce a solute distribution similar to that shown in 
Fig. 3. A fairly uniform composition specimen can 
be obtained by removing the two ends of the speci- 
men—one which is depleted in solute and the other 
enriched in solute. 

b) Reverse pass zone melting: Re-zone melting in 
the reverse direction after the initial pass will im- 
prove the uniformity of solute distribution. The ef- 
fects of reverse passing are schematically illustrated 
in Fig. 5. 

c) Use of a starting charge: By building up the 
solute concentration in the initial region of the bar, 
the initial liquid zone which forms in this region 
can be made to have its maximum solute content and 
will therefore not pick up further impurity as it is 
passed along the bar. The optimum starting concen- 
tration will be C,/k corresponding to the composition 
C, in Fig. 1. Since the liquid zone will then have its 
maximum capacity of solute, the result will be a 
composition approximating C» along the entire length 
of a metal bar, with the exception of the very last 
portion to freeze. 


Application of Zone Refining 

This paper has emphasized the basic principles 
involved in this relatively new and fascinating process 
for producing “super-pure” metals. The experimental 
techniques are relatively simple for the common 
metals with moderate melting temperatures. Zone- 
melting has been successfully applied to germanium 
and silicon for semiconductors. Work has been car- 
ried out on lead, tin, antimony, and other metals 
with astounding success. The techniques involve a 
protective atmosphere for the metal and a suitable 
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boat or crucible which will not introduce contamina- 
tion; for example, graphite boats are frequently used 
for germanium, lead, and tin refining. 

Remarkable success has been reported in removing 
arsenic from antimony by zone refining. After ten 
successive molten zone passes, the only impurities 
detectable by mass spectroscopy were zinc and arsenic 
in the order of one part in ten million. 

Some modifications of zone refining have been 
attempted for purifying reactive metals with high 
melting temperature, but this work is still in the early 
experimental stages. An example is the “cage” zone 
refining technique® which has been applied with some 
success to titanium. Molten titanium reacts viciously 
with common refractory materials thus preventing 
the use of conventional zone-melting techniques for 
its purification. In cage-zone refining, a suspended 
solid titanium specimen is moved vertically through 
an induction heating coil without contacting refrac- 
tory. If the cross section of the titanium bar is rec- 
tangular, then only the interior of the bar over the 
short length within the induction field may become 
molten; the corners of the specimen remain solid, 
owing to rapid local heat dissipation, to form a cage 
that confines the liquid. This cage effect is shown 
in Fig. 6. Surface tension effects and the levitating 
effect of the electromagnetic field, which depends on 
the design of the induction coil, prevent the liquid 
from spilling out.?/ Movement of the molten zone 
results in the purification effect. The fins or corners 


of the bar are not purified since they remain solid, 
but can be cut away leaving the purified core. Single 
crystals of titanium have not been prepared by this 
technique, but recent purification results appear very 
promising.® 

In general, zone-melting techniques provide pow- 
erful tools for metallurgical research. Of particular 
interest in the future will be the development of zone 
melting to purify many of the higher melting tem- 
perature metals. Also, expanded research will un- 
doubtedly be conducted on zone-refined metals to 
determine the true properties of pure metals and the 
relationship between minor impurities and structural 
imperfections. 
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Recognizing and Developing Creative Potential 


(Continued from page 1) 


be developed through experience and exercise. It has been 
said that nothing breeds success like success, and this is prob- 
ably true, but the corollary that failure breeds failure need 
not be true. If, through continued application, failures can 
be corrected, high orders of self-confidence can be developed. 
Actually, the fear of making a mistake is a very devastating 
emotional block to creative activity. People should realize 
that progress is made through failure as well as through suc- 
cess. I have had better success in training creative designers 
by helping them develop this spirit of self-confidence than 
I have in imbuing them with special design techniques or 
tricks of the trade. 

The last emotional attribute I will mention here is drive. 
Many people have indicated that they feel that this is the 
prime requisite of all creative workers.... 


Developing Creative Potential: The steps in the creative 
process have been listed as: question, observe, associate, and 
predict. And that last term, predict, covers a multitude of 
concepts, for it includes not only communicating your ideas 
to someone else, but communicating them in a tangible form 
which requires initiative, daring, self-confidence, and tre- 
mendous drive. Many of the best ideas had to be sold under 
extremely unfavorable conditions, and this is no easy job. 
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To do it requires faith, not only in yourself, but also in the 
worthwhileness of your project. Keen analysis, keen judg- 
ment, help provide and bolster this faith, and so they are as 
essential to this phase as they are to the idea-getting 
phase. ... 


Helping Others Develop Creative Potential: What can 
you do to help others realize their creative potential and ex- 
ercise it for the good of man? The management of creative 
personnel is a very complex problem, and, unfortunately, 
little concrete evidence is available to guide us. I do be- 
lieve, however, that management has five responsibilities to 
the creative person. It must stimulate, encourage, assist, rec- 
ognize and reward the creative person in order that he be 
most productive, though it is difficult to lay down any gen- 
eral rules for doing any one of these things and having them 
equally effective on all people. People vary considerably, and 
what stimulates one person may over-stimulate the next.... 

Some people can work to stringent deadlines; other people 
narrow their fields to such a degree that only commonplace 
results are obtained. One general rule might be made, how- 
ever, and that is to try to establish a permissive atmosphere, 
one that establishes psychological safety and psychological 
freedom for the workers in that area. 
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An Aid in the Preliminary Design of Prestressed Concrete Beams 


Jack R. CLANTON 
Associate Professor of Civil Engineering 


Introduction 


In a prestressed concrete 
beam it is desirable to have 
the full allowed unit com- 
pressive stress available to 
resist the bending moment 
from the superimposed load 
and at the same time to 
have little or no tensile 
stress on the beam section. 
In order to fulfill this ideal 
condition, the combination 
of the effect of prestress and 
the effect of the beam weight should give zero stress 
(or a comparatively small tensile stress) at the top 
of the beam and full allowed unit compression stress 
at the bottom. 


J. R. Clanton 


If the resultant of the initial prestress force (F;) 
is at the lower kern limit, its effect is to give zero 
F; 
A I 


stress at the top of the section, that is, 


2 
=0, from which wn, as shown in Fig. 1. The 


value of F;, as given in Fig. 1, to make the unit stress 
at the bottom of the beam equal the full allowed 
value, is obtained from the allowed unit concrete 

stress at time of prestress: ear Vad AP 
r*/y, is substituted for e. The effect of bending 
moment (Mg) due to the beam weight is counteracted 
by lowering the resultant of the tensioning force an 
amount ¢=Mg¢/F;, as shown in Fig. 2. Then the 
effect of initial prestress alone, and the location of 
the resultant force to give this effect, will be as shown 
in Fig. 3. It follows that the combined effect of 
initial prestress and bending due to beam weight 
will be as shown in Fig. 4. 

Subsequent to the prestressing, there will be a 
decrease in F; because of shortening of the beam due 
to concrete shrinkage, creep, and, in the case of 
pretensioned beams, elastic deformation. Creep of 
the highly stressed steel may also decrease the value 
of F; so that the final prestress force (F) may be 80 
to 85 per cent of F;; therefore, a good estimate of the 


, When 
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required section modulus may be made by multi- 
plying the superimposed-load bending moment by a 
factor of 1.18 to 1.25 and dividing by f.:. A section 
having a modulus at least as large as the computed 
value may then be selected from the chart, Figs. 5a 
and b. Finally, the selected section should be checked 
by computing the maximum stress intensities. 


The procedure discussed above is for the pre- 
liminary design at the section of maximum moment. 
For a complete design, other sections may need in- 
vestigation, and diagonal tension and anchorage 
conditions must be checked. In some cases, it may 
be more economical to revise the section to one that 
is unsymmetrical. In all cases of post-tensioned 
beams, the cables must be draped to give an ec- 
centricity of no more than e¢ at the end, if tensile 
stresses are not to exist at this point. In the case of 
pretensioned beams, if the effect of beam weight is 
to be canceled, there must be a variable beam depth, 
or a camber in the beam, so that the eccentricity 
can vary from e at the point of maximum moment 
to e; or less at the ends. 


The proportions of the beam sections, as given in 
Figs. 5a and b, follow in general the reeommenda- 
tions made by Guyon! as follows: 


Web thickness: h<30 in.: b’=d/7 to d/8 
h> 30 in.: b’=4 in.+d/40 


Flange thickness ¢: approximately 4/8 and not less 
than 0.1h for h of 5 ft to 8 ft. 


Width of flange b: not to exceed 0.8h and not less 

than 0.4h. 

Three-inch fillets were used in all sections for ease 
in preparing the chart, although in final design it 
will usually be desirable to detail them to a smaller 
slope carried out to the full flange width. In most of 
the sections a change in the fillet will not appreciably 
change the section properties; nevertheless, for a 
final design the computations should be based on the 
new section properties. 


1Y. Guyon, Prestressed Concrete, New York, John Wiley 
& Sons, Inc., 1953. 
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Preliminary Design Example 

Select a beam section to carry a 50-k load applied 
at the center of a 40-ft simple span. 

Ff, =5,000 psi; f-=0.4 2,000 psi; 

fec=4,000 psi; =0.55 ft; =2,200 psi; 


50 (40) 


M,= =500 k-ft. 


Assuming 15 per cent loss of prestress: 
M(1.18)12 — 500(1.18)12 
2.2 

From the appropriate part of the chart, the 
lightest section with the required flexural resistance 
has the following properties: 

h=38 in., b=18 in., b’=5 in., t=5 in., 

A=338 sq in., Z=3,240; 
We = (338/144) (0.15) =0.352 k/ft ; 


Required Z= = 3,220 in.* 


F=0.85 F;=316k; 


_ Me _ 61,569 


=9.59+4 2.27 = 11.86 in. 


70.4(12) 
372 


Check stress condition at time of prestress: 


372 372(11.86) 70.4(12) 
3,240 3,240 


= —1.10+1.36—0.26=0 
—1.10—1.36+0.26= —2.20 k/sq in. 


Check stress condition at time of superimposing 
loads: 


fx —316 316(11.86) 500(12) 


338 3,240 3,240 
= —0.944 1.16 —0.26—1.85 = —1.89 k/sq in. 
fo= —0.94—1.16+0.26+ 1.85 = +0.01 k/sq in. 


—0.26— 


The selected section is satisfactory, in that the 


= 70.4 k-ft; compression at the top of the beam is below the 
8 allowed value of 2.0 k/sq in. and the small amount 
Afei 338(2.2) of tension at the bottom of the beam is hardly of 
F = = 372 
2 significance. 
| 
& 
+ = 
e 
=> y, L e t 
fei 
Fic. 1. Prestress Force Fic. 2. Counteracting Fic. 3. Location of 
at Kern Limit Beam-Weight Moment Prestress Force 
(a) (c) 
4 
= 
EA 
<= 
foi 
Fic. 4. COMBINATION OF STRESS EFFECTS AT TIME OF PRESTRESS 
(a) Effect of prestress force. (b) Effect of beam weight. (c) Resultant effect. 
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OUTSIDE U.S. A. 


Highlights from Our Foreign Exchange Publications 


An APPROXIMATE SOLUTION OF THE PRANDTL EQUATION OF 
Fiumw Motion IN THE Bounpary Layer. P. Szymanski, 
Archiwum Mechaniki Stosowanej, VII-1, 1955. Polish. 
(English abstract-——more detailed than the following). 


A complete approximate solution of the well known 
Prandtl equation of fluid motion in the boundary layer is 
given in this paper. The proposed method of solution is 
based on that of Pohlhausen, the procedure, however, being 
considerably simplified. This eliminates tedious calculations 
and permits one to obtain a complete solution in the form of 
general equations which are easy to handle. 

Starting from the Prandtl equation we determine the ap- 
proximate velocity profile of the boundary layer by means 
of a fourth degree polynomial, the coefficients of which are 
determined from the boundary conditions and the equation 
of continuity. 

Thus, starting from the equations of the boundary layer 
and from the boundary conditions alone, we obtain the same 
expressions as those obtained by Loyciansky (Loitsianskii). 

Next, the thickness of the boundary layer is found from 
the Karman integral condition, thus permitting us to obtain 
a differential equation, which is analogous to that of Pohl- 
hausen. The integral \ (s) of this equation constitutes the 
solution of the problem. It vanishes for s=s,,, for which the 
potential velocity U(s) reaches its maximum, and has a 
finite value Ay for s=0. 

From this solution, we obtain the value s=sj, for which 
value the boundary layer separates. 

In spite of certain simplifications, the results of numerical 
calculations are in good agreement with those obtained by 
Pohlhausen by exact integration. 


On THE Heat Dissipation oF 4-StroKE CycLe ENGINES. 
Osamu Hiras, Report of the Institute of Industrial Science, 
Uni. of Tokyo, August, 1954. English. 


The author found that the amount of heat dissipated to 
the water jacket of a gasoline engine increased whenever 
operating conditions or design were changed in such manner 
as to increase the temperature of the combustion gases, the 
time of contact of the hot gas with. the cylinder walls, the 
area of the cylinder walls or the heat transfer coefficient be- 
tween the walls and the hot gas. A theoretical analysis sup- 
plied an equation which correlated with the experimental 
results. 


New EL-ectrotytic EtcHANTS FOR STAINLESS STEELS CoNn- 
TAINING SIGMA PuHaseE: I. Bertetti, Metallurgia Italiana, 
Vol. 68, No. 7, 1956, Italian. 


A series of new electrolytic etchants for differentiating 
sigma phase from the matrix in stainless steels of both the 
ferritic and austenitic types have been studied. For ferritic 
stainless, the optimum concentrations and times of etching 
for solutions at room temperature are as follows: 


a) 5% ammonium persulphate 5 seconds 
b) saturated solution of sodium perborate 1 minute 
c) 2% sodium periodate 10 seconds 
d) 5% sodium bromate 5 seconds 


All of these solutions darken and outline sigma phase. 
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For austenitic stainless, the periodate and perborate solu- 
tions used above do not give sufficient distinction between the 
sigma phase and any carbide phase present. The persulfate 
solution attacks and colors both the sigma and carbide. It is 
possible, however, to recognize the sigma both by the shape 
of the crystallites and the fact that they appear blue under 
the oil immersion objective. This coloration is not given 
to the carbide. A two-minute etching time is used. With the 
bromate solution, greater attack on the sigma phase than on 
the carbide is observed, permitting easy identification. An 
etching time of two minutes and 30 seconds is used. 

In all studies an anode-cathode distance of 8 mm was 
used. The current was 0.5A and voltage 6V. (Specimen 
size and anode material are not given)——Abstracted and 
translated from the Italian by E. C. ROBERTS, Associate 
Professor of Metallurgical Engineering. 


CONSIDER ENGINEERING TEACHING 


(Continued from page 2) 


be earned through summer employment with industry or 
on research contracts. These engagements are now the 
rule rather than the exception. Also, as one acquires a rep- 
utation in his field, his services are sought on a consulting 
basis. 

There is no question but what one must be prepared to 
accept some sacrifice in salary by selecting a teaching career 
rather than an industrial career—but the sacrifice is offset 
by the compensating factors noted in the preceding para- 
graphs. Happiness and satisfaction are somewhat intangible 
qualities, but the wise man will give them a high evalu- 
ation on his balance sheet. Not many college teachers 
have ulcers! 


DIGITAL COMPUTER FOR HYDRO-STORAGE 
PROBLEMS 
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with a certain amount of programming experience. 
To accomplish most efficiently the solution of large 
system studies, such as that of the entire Columbia 
River system, the use of a yet larger computer would 
be recommended. Even such problems, however, 
could probably be solved by use of a medium-sized 
computer facility. 

The method proposed in this paper is one which is 
generally adaptable for many different systems and 
types of exact solutions. More extensive work and 
further experience with the problems of hydro stor- 
age, coupled with the use of the IBM 650 computer, 
will no doubt permit further improvement and exten- 
sion of the application. 
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